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Abstract 

We investigate in detail the 77 fusion production mechanisms of two neutral Higgs bosons 
{h°A°, H°A°, h°H° and H°H°) within the framework of the mSUGRA-inspired minimal su- 
persymmetric standard model(MSSM) at an e + e" linear colliders, which provide a probe of the 
trilinear Higgs self-couplings. We calculate the dependence of the production rates on Higgs 
boson masses, the ratio of the vacuum expectation values tan j3 and the CMS energy y/s. We 
find that the cross section for the H°H° production at LC can reach 0.2 fb, while the cross 
section of A°H° production is only 10 -4 ~ 10 -3 fb under our parameters. 
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I Introduction 

Many efforts have been devoted to searching for Higgs bosons and the new physics beyond the 
standard model (SM) ^H], among which the minimal supersymmetric standard model (MSSM) 
[3] is the most promising one. Five physical Higgs bosons are predicted: one CP-odd Higgs boson 
(A ), two CP-even Higgs bosons (h° and H°) and two charged Higgs bosons (H + and H~). Until 
now all those Higgs bosons haven't been directly explored yet, only LEP2 group presents the 
strongest lower mass limits of 91.0 GeV and 91.9 GeV for the light CP-even and the CP-odd 
neutral Higgs bosons h and ^4°. For a top quark mass less than or equal to 174.3 GeV, LEP2 
group excludes the range of 0.5 < tan [5 < 2.4 4 with assuming the stop quark mixing is maximal 
and the conservative values for other SUSY parameters which affect the Higgs sector, 

The future linear colliders will continue the work in searching for physical Higgs bosons. There 
are already some detailed designs of linear colliders, such as NLCjSKy's = 500 GeV, integrated lumi- 
nosity 220fb- 1 yr- 1 ), JLC 6 \(^/s = 500GeU, integrated luminosity 220 fb^yr' 1 ), TESLA 7 \(y/a = 
500 GeV, integrated luminosity 340 fb^yr' 1 ) and CLIC 8 (1 TeV < ^fs < 5 TeV, with a lumi- 
nosity of 10 35 cm~ 2 s~ l at yfs = 3 TeV). An e + e~ LC can also be designed to operate as a 77 
collider. This is achieved by using Compton backscattered photons in the scattering of intense laser 
photons on the initial e + e~ beams ;9 ; . The resulting 7 — 7 center of mass system (CMS) energy is 
peaked at about 0.8-y/s for the appropriate choices of machine parameters. In 77 collision mode 
at the high energy peak, we may get approximately the same luminosity as that of e + e~ collision. 
Therefore, a photon LC provides additional opportunities for hunting Higgs bosons. 
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As we know, the trilinear Higgs boson couplings can be probed in Higgs boson pair production 
processes. Detailed examinations of Higgs boson production and decay processes are necessary in 
order to detect and distinguish the signals of Higgs bosons from background. The Z°h°, A°W+, 
H+W+ associated production processes at LC were studied in Refs.^3 ^2 HH El OS]- 

which have advantages in searching for heavy Higgs bosons. The h°A° pair production at one-loop 
level via e + e~ collisions has been discussed in Ref. |17|. The neutral Higgs boson pair productions 
e+e" -» h°h ,H H ,H h ,A A° and e+e" -> Z°h°h°,vi>h h at LC were studied by A. Djouadi 
et al, [HE]. The neutral ffi gg s boson P air productions (H°H°, H°h°, h°h°, A A ) at the LHC 
were investigated in Refs. pUll21M22| . The cross section for SM neutral Higgs boson pair production 
in 77 collision has been evaluated in Ref. |23| . The 77 — ► h°h° and 77 — ► A A pair productions in 
the THDM and MSSM has been calculated respectively in Ref. [21] and [231251. Ref. [23] concludes 
that the cross section of e + e~ -» 77 -» h°h° process in the framework of the MSSM at LC depends 
on m^o, tan/3, photon collision modes, and the mixing between stops. While in Ref. |26| it is found 
that for an usual mSUGRA set of parameters, 17(77 * A A ) ~ (0.1 — 0.2)fb is expected. 

In this paper we study the loop induced h°A°, H°A°, h°H° and H°H° pair productions via 77 
collisions at LC in the mSUGRA-inspired MSSM. We arrange the paper as follows: In section II, 
we present the analytical calculation. Numerical results and discussions are given in section III. 
Section IV is a short summary. 
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II Cross Section Calculation 

In our calculation we use the 't Hooft-Feynman gauge. In the loop diagram calculation we adopt 
the definitions of one-loop integral functions in reference [27]. The numerical calculation of the 
vector and tensor loop integral functions are broken down to scalar loop integrals (Ref.[2H]). The 
Feynman diagrams and the relevant amplitudes are generated by FeynArts package automatically 
|29| . The numerical calculation of the loop integrals are implemented by Mathematica package. 
II. 1 Calculation of the subprocess 77 — > (fiifo 

We denote the subprocess under investigation as 

7(&i) +7(^2) -+<h (Pi) + MP2), (2.1) 

where </>i(/>2=/i yl , H°A°, h°H° and H°H°, k\, &2, Pi and P2 are the momenta of the incoming 
photons and outgoing Higgs bosons. As the subprocess 77 — > (j>i<p2 1S loop-induced at the lowest 
order, the one-loop order calculation can be carried out by summing all unrenormalized reducible 
and irreducible one-loop diagrams and the results will be finite and gauge invariant. 

We give the Feynman diagrams of the HiA° and H{H production processes in Fig^ and 
Fig|5] respectively, where Hi(i = 1,2) represents h° and H°. The possible corresponding Feynman 
diagrams created by exchanging the initial photons or the final Higgs bosons, are also involved in 
our calculation. We find that the contributions from Z° exchanging s-channel Feynman diagrams 
with fermion loops are very small because of the Furry theorem, so we do not show them in Fig^ 
and remove them from our calculation. We also neglect the squark and the slepton loop diagrams in 
Fig^because of their vanishing contributions. That is because the contribution to the subprocesses 
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77 - ► h A , i?°^4° from the diagram with a squark(slepton) in loop which goes clockwise cancels 
exactly with that which goes counterclockwise. 

We denoted 9 as the scattering angle between one of the photons and one of the final Higgs 
bosons. Then in the center of mass system we express all the four-momenta of the initial and final 
particles by means of the 77 CMS energy \/I and the scattering angle 0. The four- momentum 
components (E,p x ,p y ,p z ) of final particles (pi and <p2 can be written as 

P 1 = { E h >\l E li~ m li sin fl > °> V E ^ ~ m ^ COS 6 ) ' 

P2 = (E^-y/El-mfamO, 0, - ^ - m*, cos *) , (2.2) 

where 

E 4>i = 2 + ( m $i ~ ' E <h = 2 ~ ( m *i ~ m ^ 2 )/^) • 

With above expressions for the four- momenta of the final state Higgs bosons in Ea. (|2.2[) . we can 
get that the three-momenta of the Higgs bosons satisfy the following relation 

\Pi\ 2 = \P2\ 2 = E 2 ^ - m\ x = El 2 - ml 2 = - [s + (m^ - m^f/s - 2{m\ + mjj] . 

The 4-momenta of the initial photons k\ and k 2 are 

h=(^V§,0,0,^V§\, k 2 = Q\/I, 0,0,- ^y/SJ. 

The Mandelstam variables are defined as 

s = (h + k 2 ) 2 = (pi +P2) 2 , 

i = (fci -Pi) 2 = (k 2 -P2) 2 , 
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u = (h - p 2 ) 2 = (k 2 - pi) 2 . 

With the definitions above, we obtain the amplitude of the subprocess 77 — > 4>i4> 2 with simple 
form. For the final states of h°A°(or H°A°), the explicit expression of the amplitude can be written 



as: 



+ + /{^J* + / 8 W « + + /{Jpfrf + fS^k la k 2p 

+ fV^haPw + f§e a ^k la p 2/3 + ffl^ktoPip + f$e a ^k 2aP2 p 
+ fSe^hahpp!^ + f§e a ^k la k 2(S p 2a kl 

+ /ge^fcia^pi^} , (2.3) 

where represents /i° or and = 1,...,18, i = 1 for i = 2 for are form 

factors. The expressions of these form factors for the two subprocesses are very similar to each other 
except some interaction vertices. We do not list their explicit expressions in this paper because 
they are very lengthy. 1 

For the final states of h°H° (or H°H°), the explicit expressions of the amplitudes can be written 

as: 

M*** = ^(fciMk) {/< V + A k) m + / s (fc) *?P? + ffKvl + /b^pTp? + / 6 (fe) « 

+ if } « + fi k) Kp v 2 + fWtftf + /gV^ + fS ] e a ^k la k 2p + fW^kujnp 
+ f[^e a ^k 2a p w + f^e^hak^p^ + f^^hak^PiaP^} , (2.4) 



lr ttie Mathematica program codes of all the form factors for 77 -> /i°A ,iy°A and 77 -> h H°,H°H are 
obtainable by sending email to zhouyj@mail.ustc.edu.cn . 



where fj(j = 1, ...15, k = 1 for h°H°, k = 2 for H°H°) are form factors. 
The total cross section for 77 — ► (ftifo can be expressed in the form 



a(s) 



1 
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t+ 



(2.5) 



spin 



In the above equation, i 



(m\ x +ml 2 -s)±yj (m^ + m| 2 - s) 2 - 4m 2 ^ J /2 , and the 
bar over summation means to take average over the initial polarizations. For the subprocess 
77 — ► H°H°, an additional factor ^ is multiplied due to the identical final particles. 

II. 2 Cross section of e + e~ — > 77 — > 0102 process at LC 

By integrating over the photon luminosity in an e + e~ linear collider, the total cross section of the 
process e + e~ — > 77 — ► (fiifo can be obtained in the form 

dC~, 



<J{s) 



dz- 



-77 



dz 



(7(77 — > (j)\(j)2 at s = z s) 



(2.6) 



where E m i n = + m^ 2 , and y/s(y/l) being the e + e (77) CMS energy. — ^ is the distribution 
function of photon luminosity, which is defined as: 



dC 
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2z I —F l/e (x)F l/e (z 2 /x) 

' Z J Xmax 



(2.7) 



For the initial unpolarized electrons and laser photon beams, the energy spectrum of the back 
scattered photon is given by j3U] 



1 



D(0 
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l-x £(l-s) £ 2 (l-x) 
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where 



D{i) = {l----,)\n(l+i) + - + - 



2 £ 2(1+0 



2 • 



(2.9) 



e = ^°. (2.10) 

m e and Eq are the mass and energy of the electron respectively, uq is the laser-photon energy, and 
x represents the fraction of the energy of the incident electron carried by the backscattered photon. 
In our evaluation, we choose ojq such that it maximizes the backscattered photon energy without 
spoiling the luminosity through e + e~ pair creation. Then we have £ = 2(1 + \/2), x max ~ 0.83, 
and D(£) = 1.8397, as used in Ref.^H- 

III Numerical results and discussions 
III.l Input parameters 

In the following numerical evaluation, we present the results of the total cross sections of the neutral 
Higgs boson pair production at LC. The SM parameters are taken as : m u = 66 MeV, = 66 
MeV, m c = 1.35 GeV, m s = 0.15 GeV, m t = 174.3 GeV, m b = 4.3 GeV, m z = 91.1876 GeV, 
m w = 80.423 GeV and a = 1/127.9|B5|. 

We take the MSSM parameters constrained in the framework of the mSUGRA-inspired scenario 
as implemented in the program package SuSpect2.1|33|. In this scenario, only five sypersymmetry 
parameters should be input, namely Mq, Mi/%, Aq, sign of \i and tan/3, where Mo, M1/2 and are 
the scalar mass at GUT scale, the universal gaugino mass and the trilinear soft breaking parameter 
in the superpotential terms, respectively. In the package SuSpect2.1 all the soft SUSY breaking 
parameters at the weak scale are obtained through Renormalization Group Equations(RGE's) |34j . 
It uses 2-loop RGE's for the gauge, Yukawa and gaugino couplings and all other RGE's are at 
1-loop order. In this work, we take Mi/ 2 =130 GeV, Ao=200 GeV and \i > 0. Mq is obtained 
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quantitatively from the m^o(or rriffo) and tan/3 values. 
III. 2 Discussion and analysis 

We depict the curves of the production rates of the processes e + e~ — ► 77 — > </>i02 i n Fig l5ITTl We 
note that if we take m^o = 200 GeV (01 m H o = 200 GeV and 300 GeV) with our chosen input 
mSUGRA parameters(M 1/2 = 130 GeV, A = 200 GeV and p > 0), the mSUGRA-inspired MSSM 
does not permit tan/3 to be smaller than 3 due to the experimental lower bound on the lightest 
Higgs boson mass m^o in MSSM (m^o > 91.0 GeV) [I]. Therefore, we take tan/? varying in the 
range of 4 < tan/3 < 38 with the above input parameters. For the same reason, we also exclude 
other parameter space regions which are inconsistent with the experimental lower bound on the 
m^o. That is why in FigEJand Fig|S]the mass of Higgs boson ^4° varies from the value of 190 GeV, 
and the curves in Fig|7]and Fig|§]start from the Higgs boson mass rrijjo = 195 GeV(for tan/3 = 4) 
and m H o = 200 GeV (for tan/3 = 30) separately. 

FigEBElare the figures for the processes e + e~ —* 77 — > h°A°, H°A° at a LC. FigEl displays the 
cross sections of the process e + e~ — ► 77 — ► h°A° as the functions of the CP-odd Higgs boson ^4° 
mass (or m h o) with tan/3 = 4 and 30 respectively with e + e~ colliding energy yfs = 1 TeV. The 
figure shows there are sophisticated structures on both two curves, which are mainly induced by 
the resonance effect. For example, at the vicinities of m^o = m ++m + ~ 300 GeV{iov tan /? = 4), 
mAo = m + + m + ~ 290 GeV(ior tan/3 = 30), m^o = 2m + ~ 570 GeV(for tan/3 = 4), and 
rrijifi = 2 m x + ~ 430 GeV (for tan/3 = 30), we can see there are visible peaks on the two curves due 
to the resonance effect from loop diagrams. 
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We plot the total cross sections of h°A° pair productions at an e + e~ LC with y/s = 1 TeV 
as the functions of tan/3 (or m h o) in Fig0J There A mass is taken to be m^a = 200 GeV and 
300 GeV respectively. The curve for m^o = 200 GeV shows that in the region of tan f3 < 5, the 
cross section goes down rapidly with tan (3 increasing, and then goes up with tan (3 varying from 
5 to 38. While the curve for m^o = 300 GeV shows the cross section decreases in the range of 
tan/3 < 5, but goes up smoothly when tan/3 > 5. That is because the coupling strengths between 
Higgs boson and quark pair are related with the ratio of the vacuum expectation values. 

Fig |21 displays the cross sections of the process e + e~ — > 77 — > H°A° as the functions of the A 
mass with tan /3 = 4 and 30, respectively. The figure shows that the cross section decreases with 
the increment of the mass of the CP-odd Higgs boson ^4°. We can see that the large tan/3 (i.e. 
tan /3 = 30) can enhance the cross section, while the moderate value of tan (3 may suppress the 
cross section. 

FigEl demonstrates the cross sections of the H°A° pair production process e + e~ — » 77 — > 
at a LC with the CMS energy yfs = 1 TeV , as the functions of tan (3 with m^o = 200 and 300 GeV, 
respectively. In the value range of tan (3 < 5, both curves for m A o = 200 GeV and 300 GeV decrease 
quickly, while in the range of tan (3 > 5 they go up with the increment of tan (3. In the region of 
tan/3 > 5, the curve for m^o = 200 GeV increases more quickly than that for m^o = 300 GeV 
with the increment of tan (3. 

Figi mTm are the figures for the processes e + e~ — > 77 — > h°H°, H°H° at a LC. Fig[7] shows 
the relationship between the cross section of the parent process e + e~ — > 77 — > h°H° and mjyo 
with y/s = 1 TeV. The large suppression on the curve of tan/3 = 4 around the position of 
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rriffo = m + + m + ~ 300GeV , corresponds to the resonance effect. Again due to the resonance 
effect, the curve of tan/3 = 30 shows a large suppression at m^o = mf x + rrif 2 ~ 205 GeV and a 
large enhancement at m^o = 2mt ~ 350 GeV. We can see from the figure that the cross section 
with tan /3 = 4 is larger than the corresponding one with tan (3 = 30 in most of the parameter 
region. 

In Fig|HJ the curves show that the cross sections of the process e + e~ — > 77 — ► h°H° decrease 
with the increment of tan (3 in the region of tan (3 < 10 (for m#o = 200 GeV) and tan (3 < 15 (for 
m#o = 300 GeV), and then their values become to be not sensitive to tan/3. 

Fig|U] indicates the dependence of the cross section of H°H° production process at LC on the 
mass of H° with Vs = 1 TeV. We find that the cross section decreases quickly with the increment 
of the mass of H°. Fig llOl shows the dependence of the cross section of e+e" -> 77 -» H°H° on 
tan (3 with m^o = 200 GeV and 300 GeV respectively. Both curves show that the cross sections go 
down at the low tan [3 region with the increment of tan /3, and it changes smoothly in the moderate 
tan/3 range. When tan/3 > 25, the curve for m^o = 300 GeV goes up slowly while the curve for 
m^o = 200 GeV keeps unchanged. The cross section of H°H° pair production at e + e~ LC is in 
the order of 10 _1 to 10 -2 in our chosen parameter space, and specially when m^jo = 195 GeV and 
tan/3 = 4, the cross section reaches 0.13/6. 

We also plot the total cross sections of the h°A°, H°A°, h°H° and H°H° production processes 
as the functions of yfs at LC in Fie llll There we take the input parameters M1/2 = 130 GeV , 
Aq = 200 GeV, fi > 0, tan [3 = 4 and m^o = 200 GeV . These input parameters induce m/jO = 
92.31 GeV and m A o = 193.6 GeV in the mSUGRA-inspired MSSM. From this figure we find 
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that the cross section for H°H° production at LC is larger than the cross sections of other three 
channels of Higgs boson pair production processes. The cross section a(e + e~ — ► 77 — > H°H°) is of 
the order 10 -2 ~ 10 _1 , and can exceeds 0.2 fb in some parameter space. It shows that this process 
could be visible at linear colliders with high enough luminosities. The curve for the cross section 
of e + e~ — ► 77 — > H°A° has the lowest value among the four curves in our chosen parameter space, 
which has the values in the order of 1CT 4 ~ 10~ 3 fb. 

IV Summary 

In this paper, We evaluate the loop induced production processes of two neutral Higgs bosons via 
77 fusion, i.e. e + e~ — > 77 — > 4>i(j)2 {4>\4>2 = h°A ,H A ,h°H and H°H°), within the framework 
of the mSUGRA-inspired minimal supersymmetric standard model at e + e~ linear collider. We 
analyze the dependence of the total cross sections on Higgs boson masses m^^i = H°,A°), the 
ratio of the vacuum expectation values tan (3 and colliding energy ^fs, respectively. With the same 
input parameters in Higgs sector(i.e. tan/3 = 4, m H o = 200 GeV), we find that the H°H°-paii 
production at LC (e + e~ — > 77 — ► H°H°) has larger production rate than the other three channels 
of the neutral Higgs boson pair productions in general, and can reach 0.2 fb quantitatively in our 
chosen parameter space, which could be visible for LC with high enough luminosity. While the cross 
sections of e + e~ — > 77 — > h°A ,H A ,h°H at LC have smaller cross sections. The H°H°-paiT 
production has the lowest production rate among all the four channels in our chosen parameter 
space, which has the values of the order 10 -4 ~ 10~ 3 fb. 
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Figure 1: The Feynman diagrams of the subprocess 77 — > Hi A , where Hi denotes h° and H°, f 
can be quark, lepton or chargino, and / is squark or slepton. The meaning of other notations in the 
figures are as follows: (a): siS2S$s±=G-G-G-H- ,G-G~w-H-, w-G~G-H-, w-G~w-H-, g-w~G-h-, 

G-W-W-H-, W-W-G-H-, W-W-W-H-, H~ H~ H~G~ , H-H-H~W-; {b):SiS2S 3 S^=G~ H~ H~G~ , 
G-H-H-W-, W-H-H-G-, W~H-H~W- , H~G~G~H~ , H-G~W-H-, H-W~G-H-, H-W-W~H-; 
(d):when the propagator is G°, siS2Ss=u±u±u±, G~G~W~, W~W~G~, and when the propagator 
is Z°, SiS 2 S3=U±U±U±, H-H-H-, G~G~G~, W~W~W~ , G-G~W-, G-W~G- ,W~G~W- ,W-W~G-\ (e): 
SiS 2 S 3 =W-G-H-, H-H-W-, W-W-H-, (f): SiS2S 3 =H-H-G-, G~G~H~, W~W-H~, H~H-W~. 
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(k) (m) 

Figure 2: The Feynman diagrams of the subprocess 77 — > HiH°, where Hi denotes h° and if , / 
can be quark, lepton or chargino, / is squark or slepton, and u denotes u+ or m_ (the ghost field of W 
boson). The figures form (a) to (j) are for both the h°H° and H°H° production processes, while (k) 
and (m) are for the H°H° final states only. si(i = 1, ...4) can be H~,G~ or W~, and they appear in the 
figures in the following way: (a):there are two cases permitted: (i)si = s 2 = S3; (ii) sis 2 S3 =G~VK~G~ , 
G-w-w-, G-G-W-, W-W-G-, W-G-G-, w-G~w-. (b):there are five cases permitted: s\ = s 2 = 

S3 = S4; (ii)Si = S2 7^ S3 = S4; (iii)Si = S4, and S2S3 =G~W~ or W~G~, (iv)S2 = S3, and S1S4 =G~W~ 
or W~G~. (v)SiS 4 =G~W- or W~G~, and S 2 S 3 =G~W~ or W~G~. (c):SiS 2 S 3 =, H~ H~ H~ , G~G~G~ , 
W-W-W-, G-G-W-, G-W-G-, W-G-W-, G-W~W- . (d): S1S2S3 =H~ H~W~ , G~G-W~, W~G~H-, 
W~G~G-, G-W-W-,W-W-H-,W~W-G-. (i): S1S2S3 = H~ H~ H~ , G~G~G~, G~W~G~ . 
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Figure 3: The cross sections a of the process e + e — > 77 —>■ h°A°, as the functions of the mass of 
Higgs boson A (starting at m A o = 190 GeV), and tan/3 is taken as 4 and 30, respectively. 
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Figure 4: The cross sections a of the process e + e~ —> 77 —> h°A°, as the functions of tan/3. 
The mass of Higgs boson A is taken as 200 Gey (starting at tan/? = 4) and 300 Ge^(starting at 
tan/3 = 3), respectively. 
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Figure 5: The cross sections a of the process e + e — > 77 —>■ H°A°, as the functions of the mass of 
Higgs boson A (starting at m A o = 190 GeV), and tan/? is taken as 4 and 30, respectively. 
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Figure 6: The cross sections a of the process e + e~ — > 77 — > H°A°, as the functions of tan/3. 
The mass of Higgs boson A is taken as 200 Gey(starting at tan/3 = 4) and 300 GeF(starting at 
tan/3 = 3), respectively. 
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Figure 7: The cross sections a of the process e + e~ —> 77 — > h°H°, as the functions of the mass 
of Higgs boson H°, and tan/? is taken as 4(starting at m H o = 195 GeV) and 30(starting at 
m H o = 200 GeV), respectively. 
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Figure 8: The cross sections a of the process e + e — ► 77 — > h°H°, as the functions of tan (3 (starting 
at tan/3 = 4). The mass of Higgs boson ^4° is taken as 200 GeV and 300 GeV, respectively. 



26 



0.14 




' 1 ' I ' 1 ' 1 - ^ 

200 250 300 350 400 

m H „ [GeV] 

Figure 9: The cross sections a of the process e + e~ — > 77 —> H°H°, as the functions of the mass of 
Higgs boson H°, and tan/3 is taken as 4(starting at m H o = 195 GeV) and 30(starting at m H o = 200 
GeV), respectively. 
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Figure 10: The cross sections a of the process e + e — ► 77 — > H°H°, as the functions of 
tan/3(starting at tan/3 = 4). The mass of Higgs boson H° is taken as 200 GeV and 300 GeV, 
respectively. 
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Figure 11: The cross sections a of the processes e + e — > 77 — ► <j)i<t>2, as the functions of the 
CMS energy y/s. 
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